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ABSTRACT 

We report high-resolution (0.07 arcsec) near-infrared polarized intensity images of the circumstellar disk around 
the star 2MASS J16042165— 2130284 obtained with HiCIAO mounted on the Subaru 8.2 m telescope. We present 
our //-band data, which clearly exhibit a resolved, face-on disk with a large inner hole for the first time at infrared 
wavelengths. We detect the centrosymmetric polarization pattern in the circumstellar material as has been observed 
in other disks. Elliptical fitting gives the semimajor axis, semiminor axis, and position angle (RA.) of the disk 
as 63 AU, 62 AU, and —14°, respectively. The disk is asymmetric, with one dip located at P.A.s of ~85°. Our 
observed disk size agrees well with a previous study of dust and CO emission at submillimeter wavelength with 
Submillimeter Array. Hence, the near-infrared light is interpreted as scattered light reflected from the inner edge of 
the disk. Our observations also detect an elongated arc (50 AU) extending over the disk inner hole. It emanates at 
the inner edge of the western side of the disk, extending inward first, then curving to the northeast. We discuss the 
possibility that the inner hole, the dip, and the arc that we have observed may be related to the existence of unseen 
bodies within the disk. 

Key words: planetary systems - protoplanetary disks - stars; pre-main sequence - techniques: polarimetric 


1. INTRODUCTION 

Planets are believed to form in protoplanetary disks. To 
study the structures of protoplanetary disks, multi-wavelength 
observations have targeted such disks around young stellar 
objects located in nearby star-forming regions such as the 


* Based on data collected at the Subaru Telescope, which is operated by the 
National Astronomical Observatory of Japan. 


Taurus, p Oph, and Upper Scorpius regions (e.g., Kenyon & 
Hartmann 1995; Preibisch & Zinnecker 1999; Furlan et al. 
2006; Andrews et al. 2009; Dahm & Carpenter 2009; McClure 
2009; Andrews et al. 2011; Mathews et al. 2012a). Among such 
nearby star-forming regions. Upper Scorpius, which is located 
at a distance of 145 pc (de Zeeuw et al. 1999), provides a unique 
environment, particularly in terms of its age. While the age is 
only ~5-10 Myr (Pecaut et al. 2012; Song et al. 2012), the lack 
of dense molecular material and embedded young stellar objects 
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indicate that the star formation process in Upper Scorpius is 
over. The area is essentially free of dense gas and dust clouds, 
and the association members show only moderate extinctions 
(Ay < 2 mag). In fact, there are only two K and M stars (1.5%), 
including the star discussed here, in the Upper Scorpius region 
which exhibit a X-band excess, although about half of the stars 
in Taurus exhibit a X-band excess (Strom et al. 1989). This 
implies that there are very few stars still hosting protoplanetary 
disks in Upper Scorpius star-forming region. 

Here, we introduce our target 2MASS J16042165— 2130284 
(hereafter J1604— 2130) in Upper Scorpius. J1604— 2130 was 
first identified as a member of the Upper Scorpius OB associ- 
ation in the spectroscopic survey of X-ray-selected sources of 
Preibisch & Zinnecker (1999), who found: stellar age 3.7 Myr, 
spectral type K2, stellar mass 1 M 0 , log 7' c n[K I = 3.658, and 
log[L/L 0 ] = —0.118. The spectral energy distribution (SED) 
of J1604— 2130 has the features of a transitional disk (Dahm & 
Carpenter 2009), in which small dust grains have been partly 
cleared from the inner disk but are still abundant at larger radii. 
Recently, Mathews et al. (2012b) presented the Submillimeter 
Array (SMA) 880 fim images, which directly resolved a face- 
on protoplanetary disk around J1604— 2130, and derived a dust 
mass of 0.1 Mj up . The availability of radio data was one advan- 
tage to selecting this target; however, other advantages abound 
and are as follows: (1) J 1604— 2130 is located in a nearby star- 
forming region, Upper Sco; (2) it hosts one of the largest inner 
holes in transitional disks discovered to date; (3) it hosts the 
most massive disk in Upper Scorpius; (4) the inclination 6° ± 1?5 
(Mathews et al. 2012b) is very close to face-on. 

Infrared photometry and spectroscopy with Spitzer has shown 
that less than ~10% of stars with disks show some degree of 
inner disk depletion yet still retain massive outer disks (e.g., 
Cieza et al. 2010; Muzerolle et al. 2010; Luhman et al. 201 1). 
Although there have been several transitional disks detected 
with radio observations with SMA around single young solar 
mass sources (e.g., Pietu et al. 2006; Andrews et al. 2011), 
only a handful of transitional disks have been detected in the 
near-infrared (e.g., Fukagawa et al. 2006; Thalmann et al. 2010; 
Hashimoto et al. 201 1). This is mainly because it is difficult to 
resolve structures with only a few 10 AU scale close to a bright 
central star. Furthermore, structures inside the hole of disks have 
rarely been directly imaged or resolved to date. To gain a better 
understanding of the morphological structures at radii where 
planets are expected to form, we conducted high-resolution 
near-infrared polarimetric observations of J 1604— 2130. 

2. OBSERVATIONS AND DATA REDUCTION 

We carried out polarimetry in H band (1.6 p m ) toward 
J1604— 2130 using the high-resolution imaging instrument 
HiCIAO (Tamura et al. 2006; Hodapp et al. 2006) with a 
dual-beam polarimeter at the Subaru 8.2 m Telescope on 2012 
April 1 1 . These observations are part of the ongoing high- 
contrast imaging survey (SEEDS; Tamura 2009). The polari- 
metric observation mode acquires o-rays and e-rays simultane- 
ously, and images a field of view of 10" x 20" with a pixel scale 
of 9.53 mas pixel -1 . J1604— 2130 was observed without an oc- 
culting mask in order to image the inner most region around the 
central star. The exposures were performed at four position an- 
gles (P.A.s) of the half-wave plate, with a sequence of PA. = 0°, 
45°, 22° 5, and 67? 5 to measure the Stokes parameters. The inte- 
gration time per wave plate position was 15 s and we obtained 25 
wave plate cycles. The adaptive optics system (A0188; Hayano 
et al. 2010) provides a diffraction limited and mostly stable stel- 


Table 1 

The Result of the Ellipse Fitting for the Disk of J1604— 2130 


Semimajor axis 

63.4 ± 1.0 (AU) 

Semiminor axis 

62.4 ± 1.0 (AU) 

Ellipse center ( u , v) a 

(-1.1 ± 1.1, 1.5 ± 1.1) (pixel) 

Position angle 

-14° ± 11° 

Inclination 

10?2 ± 10:1 


Note. a Central position (0, 0) corresponds to the stellar position. 


lar point spread function (PSF). The total integration time of 
the polarization intensity (hereafter PI) image was 750 s af- 
ter removing low-quality images with large FWHMs by careful 
inspections of the stellar PSF. 

The Image Reduction and Analysis Facility (IRAF 30 ) soft- 
ware was used for all data reduction. The polarimetric data 
reduction technique is described by Hashimoto et al. (201 1) and 
Muto et al. (2012), using the standard approach for polarimetric 
differential imaging (Hinkley et al. 2009). By subtracting two 
images of extraordinary and ordinary rays at each wave plate po- 
sition, we obtained +Q, —Q, +U, and — U images, from which 
2 Q and 2U images were made by another subtraction to elim- 
inate remaining aberration. The PI was then calculated using 
PI = \J Q} + U 2 . The instrumental polarization of HiCIAO at 
the Nasmyth platform was corrected for by following Joos et al. 
(2008). 

3. RESULTS 

The 77- band PI image of J1604— 2130 is presented in 
Figure 1(a). The emission arises from dust particles mixed 
with gas in the circumstellar structures, scattering the stellar 
light. Details of the near-infrared disk around J 1604— 2130 
have for the first time emerged at 0'. , 07. The disk has a dip 
D, located around a PA. ~ 85°, making this transitional disk 
asymmetric. Inside the bright disk, we see a region where the 
surface brightness drops by a factor of ~5 compared with the 
disk. We call this inner lower-brightness region the “hole” 
throughout the remainder of this paper. Figure 1(b) shows 
77- band polarization vectors. The disk with its centrosymmetric 
polarization pattern surrounds the central source J 1604—2130 as 
previously observed in other disks using the technique of near- 
infrared imaging polarimetry (e.g., Hashimoto et al. 201 1). 

Elliptical fitting was performed to measure several disk pa- 
rameters including the angular separations of the major and 
minor axes, the P.A., inclination, and the position of the ellipse 
center. We first measured disk radial profiles in 5° P. A. in- 
crements, and then extracted coordinates of the brightest peak 
area. Those coordinates were used to fit an ellipse. The non- 
linear least-squares Levenburg-Marquardt algorithm with the 
above five free parameters was performed and its result is dis- 
played in Table 1 . Its shape is superimposed on Figure 1(a). The 
derived semimajor and semiminor axes are called “NIR peak 
radii” in this paper. The geometric center of the disk is consis- 
tent with the position of the star within errors. The disk P. A. 
has a relatively large error because the disk is close to circular. 

Figure 2 shows the azimuthal surface brightness at r = 33, 
63, 145 AU, which are illustrated in Figure 1(c), with the 
P. A. measured from north to east. The error bars shown in 
Figure 2 represent the calculated standard deviation. The left 


30 IRAF is distributed by National Optical Astronomy Observatory, which is 
operated by the Association of Universities for Research in Astronomy, Inc., 
under cooperative agreement with the National Science Foundation. 
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Figure 1. //-band HiCIAO images of J1604— 2130. The saturated central area (radius = if. 2) is masked in black, (a) The PI image of J1604— 2130. The field of view 
(FOV) is 2'.'9 x 2'.'9. The unit of the color bar is mJy arcsec” 2 . The light blue ellipse and plus sign are the best-fit result of our elliptical disk model and the ellipse 
center, (b) H - band polarization vectors superposed on the PI image. The vector directions indicate angles of polarization. The plotted vectors are based on 7 [pixel] x 
7 [pixel] binning corresponding to the spatial resolution. The FOV is 2'.'0 x 270. The vector’s lengths are arbitrary, (c) Red (r = 145 AU), brown (r = 63 AU), and 
yellow (r = 33 AU) circles, corresponding Figure 2, superimposed on the PI image, (d) SMA 880 //m continuum map (Mathews et al. 2012b) superimposed on the 
PI image. White color contours indicate 2rr. 3 a, 6a, 9a, and 1 2rr intensity (la = 1.3 mJy beam 1 ). The ~075 x 073 beam of SMA is shown in the bottom right. 



Figure 2. Azimuthal surface brightness profile at r = 33 AU (inner hole region), 
r = 63 AU (inner disk region), and r = 145 AU (outer region) with position 
angle measured from north to east. The radius 33 AU is measured as it is the 
brightest part of the arc. 

panel of Figure 3 shows the radial surface brightness profile 
of J1604— 2130 along the minor axis. The error bars shown 
in Figure 3 represent the calculated standard deviation. In the 
north, the surface brightness along the minor axis increases as 


^3.6±o.2 f rom 40 to 62 AU, and decreases as r _4 - 7±al from 74 to 
124 AU. In the south, the surface brightness along the minor axis 
increases as r 3 0±0 - 4 from 40 to 62 AU, and decreases as r ~ 4J±0 - 2 
from 74 to 124 AU. The right panel of Figure 3 shows the radial 
profile of the surface brightness along the major axis. In the east, 
the surface brightness along the major axis increases as r- 6±0 - 2 
from 40 to 62 AU, and decreases as r ~ 4 - 7±ai from 74 to 124 AU. 
In the west, the surface brightness along the major axis increases 
as r L9±a2 from 40 to 62 AU, and decreases as r -4 ' 0±a2 from 74 
to 124 AU. The profiles along both the major and minor axes 
show a change of slope beyond 63 AU, thus reconfirming the 
NIR peak radius to be 63 AU. Along all four axes, the counts 
drop to half of the peak counts at around 50 AU from the central 
star. The width of the brightest region where counts exceeds half 
of the peak counts in the disk is around 30 AU. 

A salient feature inside the hole is an elongated arc extending 
over the disk inner hole. It begins at the inner edge of the 
western side of the disk, extending inward first, then curving 
to the northeast. The whole arc extends at least 50 AU. This 
arc is equivalent in brightness to the inner edge of the disk 
and to the dip D region. This arc consists of a real feature 
and is not an artifact because in contrast to the inner hole 
errors (shown as squares in Figure 2), the surface brightness 
difference between the arc and the other regions in the inner 
hole has a 6a significance, whereas the other directions show 
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Figure 3. Left panel: radial profile along the minor axis. NW and SE radial profile are averaged over —25° < P.A. < 5° and 155° < P.A. < 185°, respectively. This 
is because the disk position angle derived from our near-infrared image and peak velocity of the sub-mm CO emission (—5° ± 10°) are 9° different. Therefore, we 
centered our 30° wide radial profiles on the average of our near-infrared and sub-mm line emission. Right panel: radial profile along the major axis. In the same way, 
NE and SW radial profiles are averaged over 65° < P.A. < 95° and 245° < P.A. < 275°, respectively. 


small fluctuations in emission. Furthermore, the polarization 
vectors in the region of the arc mostly face the central star as 
well as the region of the bright disk, indicating that the arc is 
also illuminated by central star and is physically connected to 
inner edge of the disk. 

In the same way, the dip D is a real dip feature and is not an 
artifact because the difference in surface brightness between D 
and the rest of the disk is statistically significant and not within 
error bars. In addition, the polarization vectors in the region of 
the dip mostly face the central star. 

There appears to be another dip located around a P. A. 
~255°, which is not azimuthally symmetric with D in terms 
of location and surface brightness. However, this another dip is 
not significant outside of the uncertainties of the observations, 
thus follow-up future observations are required to confirm its 
presence. 

4. DISCUSSION 

4.1. Comparison with the SMA Imaging 

Mathews et al. (2012b) present SMA 880 pm continuum and 
CO (3-2) line observations of this source with a resolution of 
0/34 and resolve the disk (Figure 1(d)). They determine the 
surface density distribution of the disk with inner cavity whose 
radius is 72 AU based on the model fitting to the visibility at 
880 /xm. They also suggest that sparse matter still remains inside 
the cavity (20-72 AU) because of the dust emission indicated by 
mid-infrared (Carpenter et al. 2009) excesses and their model 
fitting. The NIR peak disk radius of 63 AU derived from our 
H - band image roughly agrees or may be slightly smaller than 
their measurement and can be explained by scattered light 
reflected at the inner edge of the disk. 

Based on their CO peak velocity map, Mathews et al. (2012b) 
derive the P.A.s of the major and minor axes to be 85° ± 10° and 
—5° ± 10°, respectively. Their derive disk parameters, including 
the P.A. and inclination, are in a good agreement with values 
derived from our near infrared observations. 

Mathews et al. (2012b) point out that the relatively high mass 
of the gas-dust disk, its sharp inner edge, and the lack of gas 
accretion onto the star together point to a dynamical origin for 
the inner hole. No companions have been found down to brown 
dwarf masses (Kraus et al. 2008; Ireland et al. 2011). Based on 


their analysis, they suggest that the origin of the hole can be 
attributed to the unseen presence of one or several giant planets 
within about a 40 AU. 

The width and depth of the hole depend on the planet mass, 
disk viscosity, and disk thickness (e.g., Takeuchi et al. 1996; 
Crida et al. 2006). As discussed by Dong et al. (2012a, 2012b), 
the radius of the inner hole is on the order of ~60-70 AU. There 
have been extensive studies on the process of gap opening by 
a planet embedded in a disk (e.g., Lin & Papaloizou 1986; 
Takeuchi et al. 1996; Crida et al. 2006). A rough estimate of the 
upper limit of the size of the gap that one planet can make may be 
given by the location of the lowest order Lindblad resonances, 
which, in the case of the outer resonance (the outer edge of 
the gap), is at ~1.6 times the orbital radius of the embedded 
planet (Takeuchi et al. 1996). If this is the case, the existence 
of the inner hole might indicate the existence of at least one 
planet at ~40-50 AU from the central star (see also Crida et al. 
2006), and the inner edge of the hole would be expected at 
~25-30 AU, which is inside of our saturation radius and is 
therefore unresolved in this study. 

If multiple planets reside within the hole, the hole width 
can be wider and the inner disk may be much smaller than the 
single planet case (Zhu et al. 2011). To understand the properties 
of the whole of this object, it is important to give constraints 
on the inner disk structure using detailed SED modeling (see 
Dong et al. 2012a, 2012b) and/or high-resolution observations 
by ALMA. Our observations do not show a clear inner disk 
edge, suggesting that the inner disk radius must be inside of the 
saturation radius of our observations. 

4.2. Non-axisymmetric Features 

We list a few possibilities for origin of the two non- 
axisymmetric features, the dip and the arc since it is difficult 
to draw firm conclusions at this stage. According to Figure 6(c) 
in Mathews et al. (2012b), the only structures that are not re- 
produced by the axisymmetric disk model within 2 o error are 
the southeast component; the north dip is consistent with the 
axisymmetric disk model convolved with the elliptical beam. 
We, therefore, do not elaborate on the north dip on the submil- 
limeter (sub-mm) disk in the following discussions. Mathews 
et al. (2012b) observe that the millimeter continuum emission 
has a peak at O'/ 5 southeast from the star in the 880 p m image. 
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A similar distribution is also seen in the total integrated inten- 
sity distribution of CO (3-2) emission. Our //-band image and 
Figure 3 show that the southeast of the disk is brighter than the 
east and the west, despite the fact that we observe the scattered 
light originating at the disk surface. This may be explained in a 
way that the more materials exist, the thicker the disk is, result- 
ing in a higher grazing angle of the stellar light at the scattering 
surface (e.g., Muto 2011). In our near-infrared image, the north 
of the disk is also bright. The brightness of the disk in the north, 
together with sub-mm observations suggesting material accu- 
mulation in the south, infers that the northern portion of the disk 
is the near side. This assumes that forward scattering dominates. 
If this picture of the disk orientation is correct, the disk rotation 
is counterclockwise direction. This geometry might also ex- 
plain the origins of the non-axisymmetric structures at D. D is 
at the interface between where the material accumulated (south) 
and where forward scattering is significant (north). However, it 
should be noted that since the disk is nearly face-on, signifi- 
cant forward scattering excess and precise prescriptions of the 
phase function should be required, and our discussion crucially 
depends on the details of the optical properties of the dust par- 
ticles. It should also be noted that the material concentration in 
the south is within 2 a of the current sub-mm data (Mathews 
et al. 2012b). The detection of the non-axisymmetric structures 
at sub-mm wavelengths is needed to examine the details of the 
disk geometry. 

Another possibility for the origin of the dip D is that there is 
a massive gravitating object at the disk’s midplane. It could 
be a planet (which is not necessarily the same as the one 
discussed in the previous section as the origin of the in- 
ner hole), a clump formed by hydrodynamic instability (e.g., 
gravitational instability), or something else. Such objects can 
pull the materials from the surface toward the disk’s mid- 
plane, thereby producing a shadow in the scattered light image 
(Jang-Condell 2009). This might produce features such as D. 
Such objects have sizes smaller than the SMA beam size, there- 
fore cannot be resolved with the SMA. 

For the arc structure in the west, the morphology is similar 
to spiral density waves in the disk. An arguably uniqueness is 
that this arc extends over the “disk inner hole” region. It also 
resembles to spiral density waves in terms of the pitch angle. 
Spiral waves can be produced by dynamical processes in the 
disk such as the turbulence or the planet-disk interaction (see, 
e.g., Papaloizou et al. 2007). If this feature is a spiral wave, its 
morphology indicates that the disk rotation is counterclockwise 
because spirals are in general trailing features. This rotation 
direction is consistent with the picture presented at the beginning 
of this section. Last, we note that the arc structure shows 
morphological similarity to the “planet shadows” discussed 
by Jang-Condell (2009). One possible way to distinguish the 
different scenarios is to observe the structures of the disk’s 
midplane at longer wavelength. For instance, if a density 
enhancement is observed at the location of D, this might indicate 
the existence of a massive body there. If we see an arc structure 
in the sub-mm as well, then it is likely a density wave. 

Another possible way to reveal the nature of the arc may 
be the contemporaneous SED and scattered light observations. 
Previous studies indicate that substantial variability has also 
been seen in MIR SED and has been linked to variable structures 
including arcs or rims of the disk (Dahm & Carpenter 2009; 
Wisniewski et al. 2008; Muzerolle et al. 2009; Espaillat et al. 
2011; Flaherty et al. 2012). For the J1604— 2130 case, the 
location of the arc and the rim of the disk wall is generally 
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consistent with the radial location in the disk which possibly 
contributes to the SED variability. 
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